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Abstract. The global constraints on chaotic dynamics induced by the analyticity of 
^vj ■ smooth flows are used to dispense with individual periodic orbits and derive infinite 

families of exact sum rules for several simple dynamical systems. The associated 
Fredholm determinants are of particularly simple polynomial form. The theory 
^ ' developed suggests an alternative to the conventional periodic orbit theory approach 

to determining eigenspectra of transfer operators. 

o 

(N 

l>- ! PACS numbers: 0320, 0365, 0545 



AMS classification scheme numbers: 58F20 



c 

6 . 

Dedicated to Adrian Douady for his 60th birthday 

O 

K/t \ 1. Introduction 

H . 

Low dimensional chaotic classical and quantum dynamical systems EL 0] can be 

analyzed in terms of unstable periodic orbits. The periodic orbit theory of such systems 
has been successfully applied to a wide range of physical problems |3|, [|, |j. However, 
since the phase space is tessellated into linearized neighbourhoods of periodic points, 
analyticity properties of smooth flows are not fully utilized in the conventional periodic 
orbit theory. 

In this paper we explore global constraints on chaotic dynamics induced by the 
analyticity of the flow. We propose to dispense with the periodic orbits altogether, and 
extract spectra of transfer operators from dynamics-induced relations among periodic 
orbit sums. 
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For example, evaluating trace formulas for 2-dimensional Hamiltonian flows [pj we 
have found exact periodic orbit sum rules 

^ A^l-l/A,) 2 ~ 
ieFix(/») n ' l> 

v A, + 1/Aj 

^ AJl-l/A,) 2 

where / is a volume preserving Henon map and A, is the expanding eigenvalue of its 
ith periodic point, real or complex. Such relations are remarkable insofar they require 
a high degree of correlation between the stabilities of exponentially large numbers of 
periodic orbits distributed over the entire phase space. As a matter of fact, for at 
least two of the dynamical systems studied here, the quadratic polynomial map and 
the Farey map, the information carried by the periodic points is so redundant, that 
the periodic points must satisfy infinitely many independent sum rules. 

The key idea behind the method by which the sum rules are derived is to replace 
a weighted sum over periodic orbits by an integral representation; the dynamics then 
induces recursion relations between these for cycles of different topological length. In 
posing the problem this way we draw inspiration from the work of Eremenko, Levin, 
Sodin, and Yuditskii |t], ||, |) and Hatjispiros and Vivaldi [J10| on the polynomial 
mappings, and Contucci and Knauf [O for the Farey map. For these models the 



recursions relate hnitc numbers of terms, and yield expressions ( |27[ ), (|6l| ) for the 



corresponding Fredholm determinants in terms of finite polynomials of form 

det(l -zC) =dct(l-zL) 

where £ is a weighted transfer operator and L is the corresponding hnite [£x£] 
dimensional matrix that relates £ consecutive sums. 

While this fmiteness is a very special property of the particularly simple models 
considered and cannot hold for a generic dynamical flow, the method nevertheless 
suggests that one might be able to dispense with the periodic orbit theory altogether, 
and manipulate instead the traces of dynamical transfer operators by means of 
such integral representations. This tempting possibility motivates us to explore the 
applicability of the method in some depth, the main results being the sum rules and 
closed form expressions for Fredholm determinants of weighted transfer operators. 

We start with polynomial mappings in section |3|, and then specialize to quadratic 
mappings (section Et) for which an infinity of sum rules is easily derived. Special 
examples are worked out in the appendices. In section H we generalize a result of ref. [pf 
about the spectrum of a non-polynomial transfer operatorfor quadratic mappings. In 
section |6| we show that our method can be extended to higher dimensions and apply 
it to the Henon map and the kinematic dynamo model. Finally, in order to show that 
the contour integrals are not the essence of the method, in section we derive the same 
kind of explicit polynomial Fredholm determinants for the circle-map (or spin-chain) 



thermodynamics. We discuss the applicability of the method to continuous flows in 



section 3.1: the sum rules offer a new invariant characterization of the errors due to 



approximating a flow by a polynomial Poincare section return mapping. The sum 



rules of section M use a "signed" rather than the "natural" measure; in section 8.2 we 
outline how they might be applicable to the period-doubling presentation function, 
for which the "signed" measure yields a formula for the Feigenbaum 5. We finish with 
a critical summary and relegate some details to appendices. 

2. Transfer operators and Fredholm determinants 

We start with a brief summary of the dynamical systems concepts needed in what 
follows. For motivation and background we refer the reader to literature, for example 
refs. @,g. 

A dynamical flow is usually investigated via its Poincare return map /, with the 
evolution of a density function <fi(x) given by 

(£.<f>)(y) = Jdx5(x~r 1 (y))cj ) (x)= £ <j>{x) (1) 

v- f(x)=y 

(for the time being we assume that / is a 1-dimensional map). A variety of 
applications H, n2[ , some of which we shall turn to below, requires use of generalized 
transfer operators weighted by integer powers of the stability of the orbit 

C {k) (y,x) = f(x) k - l 6(x-r l (y))- (2) 

Let f n (x) be the nth iterate of the map, and 

ra-l 
3=0 

be the linear stability evaluated at the periodic point Xi , given by the product over the 
n periodic points belonging to a given n-cycle. We shall denote the nth iterate raised 
to mth power by f n (x) m , and the derivative of the nth iterate by f n '(x) — j^f n {x). 
Traces of the powers of this transfer operator are given by 

tr£ «= E tt^ttaT (4) 

where we have assumed that all cycles have stability eigenvalues A,; ^ 1 strictly 
bounded away from unity. A cycle is called attracting, neutral or repelling if |Aj| < 1, 
| Ai| = 1, | Ai| > 1, respectively If the map is repelling and all |AJ > 1, we can drop 
the absolute value brackets in (Q) and define the trace by 

Tn(k):= £ -^— . (5) 

x i6 Fix(/") v l ' 



A corresponding "determinant" or "Zeta function" is related to the traces by 



F{z, k) = det(l - zC {k) ) = exp I - J] — T n {k) = }_, a n (k)z n , (6) 

\ n=l / n=0 

and the spectrum of the transfer operator is given by the zeros of F(z,k). In 
applications [|14| Q determinants are often preferable to the trace sums. For a lack 
of a better term we shall refer to this function as the "Fredholm determinant" of the 
transfer operator (||) even though, strictly speaking, the term should be used only for 
determinants of operators proven to be trace class. 

The main difference between the generalized transfer operators defined by (0), and 
the operators that usually arise in the dynamical systems theory is that instead of the 
"natural measure" l/|A p | we deal here with a "signed" and possibly complex cycle 
weight 1/Ap. There are several important applications of transfer operator formalism 
based on such "signed" measures, in particular the periodic orbit formulas for the 



spectrum of the period-doubling operator, section 8.2. 

From now on we define the trace of an transfer operator C by the above trace sum 
(H) over all cycles, real and complex. This choice of the cycle weight is motivated by 
the contour integral sums to which we turn now. 

3. Contour integral evaluation of periodic orbits sums 

A sum that runs over zeros of a function can be cast as a contour integral residue 
calculation. Periodic orbit formulas are sums over zeros of the periodicity condition 
X — f n (x) = 0, therefore we consider contour integrals of type 

T = I — h ^ 
J 2iri f(x) - x 

where / and h are polynomials, and the contour encloses all zeros of f(x) — x. By 
shrinking the contour to a sum of contours around individual zeros, figure [l| and using 
the Cauchy formula 

dx 1 

^— : -r = Ok,i 
Zm x K 



we obtain a sum over fixed points Xi of the map / 

hjxj) 



T= E TT^V (7) 



x 4 eFix(/) 

Note that this cycle sum includes all periodic points, real as well as complex. For /", 
n — > oo these points fill out the Julia set, the closure of the set of the repelling periodic 
points |lq| . We assume that all zeros are simple; multiple zeros can be treated as well. 
Here we shall compute 




Figure 1. Shrinking the contour C to a sum of contours around individual zeros 
leads to a sum over fixed points Xi of the map /. For a sufficiently large contour 
C, the iterate C — > /(C) wraps twice around the periodic point set, leading to the 
key formula (E4h. In this figure the map is f(x) = x 2 — c, c = 0.5 — 0.7i, and the 8 
periodic points plotted are the zeros of f 3 (x) — x = 0. 



for k = 0, 1,2, 3, . . . (we turn to the k = —1 case in section pi). 

The key observation is that T can be evaluated by pushing the contour to \x\ — > oo, 
where the individual fixed points play no role. When k > 0, the following expansion 
of (@) is convenient: 



T(k) = & 7— 



dx_ j'(x) k x/fix) 
2iri 



x 1 

dx f{x) k ( 
2iri x 



x/f(x) 



m \f(x) 



m 



(9) 



For a polynomial this is an exact formula; we use the fact that for sufficiently large 
\x\ and order N > 2 of the polynomial /, (1 — x/ f(x))~ 1 can be expanded as a 
geometric series. Terms with powers of (xj f(x)) r higher than r > k are of order at 
least 0(x~ N ) and do not contribute to the \x\ — > oo contours. Conversely, if / is an 
analytic function, no such finite series representation is immediately apparent. 



By rescaling and translation any polynomial of order TV can be brought to form 
fix) — x N — cx N ~ 2 — . . . . Hence the last term is given by 

dx 1 f'(x) k _ I dx 1 N k x k( - N ~ 1 \l + 0(x- 2 )) _ . 



2m x (f(x)/x) k J 2m x x k ( N - l )(l + 0(x- 2 )) 
The first term in (g) is a sum over the zeros of f(x) 

^-/sf-J:*' (11) 

JK ' /(x)=0 

By the chain rule (0), the nth iterate sum is given by 

c n (k)= e r'(*) k - l = E f^'iy)^ 1 E /'(s)*- 1 -^) 

/»(x)=o f n -Hv)=° f(*)=y 

For k < 3, the sum over zeros of f(x) — can be replaced [J7J by the pre-images 
f{x) = y of arbitrary constant y, as for fc = 2, 3 the difference 

dx f /'(x) fc /'(*)*' 



2tti [ /(») - y /(a;) 

_ y : ,--i / ^ ( NxiN ~ 1} (£ 2 ) caJV ~ 3 - : : ■)* (13) 

f^ v J 2ni ix N - cx N - 2 - . . .y 

vanishes, and hence by induction the nth iterate sum ([12) is given by M 

C„(2) = C(2)" C„(3) = C(3) n . (14) 

The value of the sum C(k) depends on the parametrization of the particular 
polynomial. For example, for the cubic polynomial 

fix) = x 3 -cx~d (15) 

we have 

dx_ (3*»-c)» /fr^ Ajc 

2m x 6 —cx — d J 2m V 3x z x 



cm=t^^—_ b^a 32 M i -a+^+--)-3 C . 



If the polynomial is even, fix) — y(x 2 ), the terms for (x/ fix))-* for which fc + J is 
odd vanish, as 

dx f'jxf f x V : _ fe r dxx^- 1 g'jx 2 ) k 
2ri~^\Jix)J ~ J 2™ ,g(x 2 )J ' ( ' 

For example, for even polynomials C(2) = 0. For general k there are no formulas of 
the above type for the other terms in the expansion (S). 

3.1. Sum rules for periodic orbits of any polynomial 

Our periodic orbits sum rules are obtained by pushing the contour to \x\ — » oo in (ph. 
The first periodic orbits sum rule for 1-d maps follows from 

f dx 1 / dx 1 

( ' = J 2ri fix)-x = J 2m" x N il + Oix- 1 )) = ' 



hence 

£ a^T = °- (17) 

teFix(/») l 
The second and third periodic orbits sum rules follow from (101): 

7 27Ti r"(x) — x z -^ Aj — 1 

17 ^ v ' iePix(/») 

r w (2)= / j^ /" ,( ^ )2 = V -^- = TV 2 " + C(2)". (19) 

v ; / 2?ri f n (x)-x *-*> Ai-1 w v ; 

The corresponding Fredholm determinants follow from (1q) : 

F(«,0) = 1, F(z, 1) = 1 - zN, 

F(z,2) = (l-zN 2 )(l-zC(2)). (20) 

The T„(0), T"„(l) sum rules will be generalized to d-dimensions in section g. Using (g) 
we obtain the fourth sum rule 

r. m . cm- + / * ^|! + ** . j^ _£_ . (21) 

For even polynomials the second term vanishes. 

These are exceedingly simple sum rules and expressions for Fredholm determinants, 
and we can already discern both their utility and the ways in which they might fall 
short of what we need in order to perform efficient dynamical systems computations. 
By taking into account all periodic points, not just the real ones, these sum rules 
dispense with the infinite intricacies of controlling the parameter dependence of 
symbolic dynamics of systems as simple as a parabola. However, there are clearly very 
important physical effects that our sum rules are blind to: the most striking is the 
role played by bifurcations. While at a bifurcation real dynamics changes qualitatively, 
going from hyperbolic through intermittent to a stable attractor, the corresponding 
complexified "signed" sum rules are arranged precisely in such a way that these effects 
cancel. 

Starting with T(3), the sum rules cannot be cast in form applicable to arbitrary 
polynomials, so we now specialize to the quadratic polynomials. We refer the reader 
interested in the transfer operators weighted by more general rational polynomials to 
ref. Hi. 



4. Sum rules for the quadratic map 

The quadratic polynomial 

f(x) =x 2 -c (22) 



is the simplest example of a non-linear 1-d mapping (note the unconventional sign of 
the parameter c, chosen to simplify our final formulas). By the chain rule (J3|) for the 
derivative of an iterated function, the stability of an n-cycle {xi, X2, • • • ,X n } is given 

by 

A = 2 n Xl x 2 ---x n . (23) 

We shall now derive an explicit formula for the sum rules and Fredholm 
determinants for transfer operators (0) weighted by integer powers of the unstable 
eigenvalue A fe in (||) . The derivation is based on the following simple observation which 
relates contour integrals of functions of the nth iterate f n (x) to contour integrals of 
functions of f n+1 (x): for a quadratic polynomial and a sufficiently large contour, a 
change of the integration variable x — > y — f(x) wraps the contour around twice 

i d ^Km) = 2i^h{y). (24) 

The factor 2 arises because for large |x| the phase of f(x) « x 2 advances at twice the 
speed of the phase of x, see figure H. Reconsider now the sum (||) with the last term 
evaluated as in (|1C)|): 

The individual terms in the sum are the diagonal terms of the [(k— 1) X (k— 1)] matrix 

A^:=I^^^f n '(x) k , r,a=l,2,.-.,*-l ) (25) 

J 2iri f n {x) r 

so the nth iterate trace can be expressed as 

T„(fc) = trA (n) +2 fc ". 

Trivially, for even polynomials A rs 7^ if k and s—\ are either both odd or both even. 
The motivation for introducing the matrix A^ n > is the observation that its elements 
are a convenient basis for relating the successive trace sums T n (k) — > T n+ \(k). By 

© 

A (n+1) _ I dxf'(x) X s f'(x) fnlfff \\k 

vanishes for even polynomial / unless s + k — even. Hence we can substitute 
x s - 1 f(x) k - 1 = 2 k - 1 x s+k - 2 = 2 k - 1 (y + c) 5 ^ 2 , and applying (§|) obtain 

We note that this step relies on the very simple form of the quadratic polynomial, 
so a generalization to arbitrary polynomial mappings is not immediate. Now expand 
binomially the (y + c) term 

2 / fc+s-2 \ 



p=0 v y 



and observe that the successive A^ are obtained by multiplication by the [(k — 1) x 
(k — 1)] dimensional matrix 

if k, s differ in parity 

L(A)r« := { / k+s-2 \ ,, (26) 



2 

r- 1 



fc + s 

c 



otherwise. 



As A' ) = 1, the nth level trace sums are given by the trace of A^ n ) = 2 kn h? k -., and 
the Fredholm determinant (Q) of the transfer operator (g) is given by 

F(z/2 k , k) = (1 - z) dot (1 - zL (k) ) (27) 

The first 10 such determinants are tabulated in table [l| and an example is given in 



Appendix A. 5. This is the main result of this section: we have obtained an explicit 
formula for the spectrum of the transfer operator C^k) for any positive integer k in 
terms of the eigenvalues of a finite matrix L^j . As by the symmetry of / the entries in 
the half of the columns of (E6|) vanish, in computations it is convenient to distinguish 
the odd and even cases, and take 'L(k) to be a [/ x /] matrix, i, j = 0, 1, 2, • • • , I — 1 : 

l +Ac l+j ~ 2i odd k = 2l + l 

L.,,,,^-5 . 2 \ (28) 

+ J + 1 V+^ 2 * even fc = 2? + 2. 

These finite matrices were first introduced by Levin, Sodin, and Yuditskii, section 4 
of ref. @. 

Table 1. Fredholm determinants (127]) for the quadratic map transfer operators 
£(k), for k a positive integer. 



_k dot (1 - zL (fc) j 

2 1 

3 1 — cz 

4 1 - 2 c z 

5 1 - (3 + c) c z + 2 c 3 z 2 

6 1 - (4 + 3 c)cz + 8 c 3 2 2 

7 1- (5 + c)(l + c)cz+ (20 + 5c + 3c 2 ) c 3 2 2 - 8 c 6 ; 

8 1 - (6 + 10c + 4c 2 ) cz+ (40 + 24c + 20c 2 ) c 3 z 2 - 64c 6 r 

9 1- (7 + 15c+10c 2 +c 3 ) cz+ (70 + 70c + 82c 2 + 14c 3 +4C 

- (280 + 64 c + 28 c 2 + 20 c 3 ) c 6 z 3 + 64 c 10 2 
10 1- (8 + 21c + 20c 2 + 5c 3 ) cz+ (ll2 + 160 c + 250c 2 + 98 c 3 + 40 c 4 ) c 3 2 

- (896 + 512 c + 320 c 2 + 280 c 3 ) c 6 2 3 + 1024 c 10 2 



6 ^3 

4~\ ^.3 .,2 

10 ,4 
3 „2 

4 
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4-.1. Dynamical zeta functions 

Using a different approach, Hatjispiros and Vivaldi fllOfl have introduced a family 
of dynamical zeta functions for complex quadratic polynomials, and conjectured 
that these zeta functions are of a particularly simple rational form. We note that 
the formula (U% proves the conjectured rationality of this family of dynamical zeta 
functions. However, their method of evaluating the determinants is still of considerable 
interest, as it appears to be more efficient for evaluation of det (l — zLn,)) for high k 
than the direct evaluation undertaken here. 

For purposes of comparison with ref . ]10[ we need to relate dynamical zeta functions 
to the Fredholm determinants given above. For |A| > 1, we expand the weight in (|5|) 

1 1 ^ 






A r (l-1/A r ) A r 

and obtain the product representation of the Fredholm determinant for 1-d maps 

( °° z npT 1 \ °° ( z rlp \ 

P „,*,.„^j : g_ 3? ^- IJ j-nn('- S p i TJ • 

where the p product goes over all prime (not self-retracing) cycles. The dynamical 
zeta function is defined Q as 

vc (fc ) = n ( i - ^i) > (M) 



p 



H< 



so for 1-d mappings it can be expressed as a ratio of two Fredholm determinants 

l/C ik) =F(z,k)/F(z,k + l). (30) 

The zeta functions of ref. JlO] are indexed and normalized differently; here we follow 
the conventions used in ref. |1^] . Our Fredholm determinants are related to the 
polynomials G of ref. |1(| by det (l — zL/m) = Gk—i(z/2, — c), and our explicit formula 
(p?j) is a proof of their conjectures. 

5. Julia set escape rate 

So far we have obtained explicit finite expressions for the spectral determinants of 
signed transfer operators C(k) with weight (f') k ,k > 0. These results can easily 
be generalized to arbitrary polynomial weights . However, the methods outlined in 
section || fail for weights (f') k with k < 0. Especially the case k = — 1 is physically 
interesting, because the largest eigenvalue of the operator gives the escape rate 7 from 
a given enclosure T, which is defined by the fraction of initial points that stay in T 
after n iterations Q 

_ n ^ J r dxdy6(x-P(y)) 
6 ~ J r d X ■ (31) 
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The signed version of this transfer operator is defined with the weight ,,} , 2 and by 
summing over all real and complex pre-images of x: 



(C ( ^)(x) = £ j^L, x,yeC. (32) 

If the Julia set is real the two operators coincide. 

The spectral determinant of J32| ) has been computed in ref. g for quadratic maps 
f{x) = x 2 — c. If c lies outside the Mandelbrot set M, the critical point is attracted 
by infinity and the spectral determinant is given by the entire function 



F(z ' 1)_1+ §2«/(o)/ 2 (o)...no)' (33) 

which involves only forward iterates of the critical point of /. If c lies inside the 
Mandelbrot set, the critical point is attracted by a limit cycle. If the stability of this 
cycle is different from zero, (|33| ) is the form of the spectral determinant for small z and 
can be meromorphically continued. If on the other hand the critical point is periodic 
with period n p , the spectral determinant is a polynomial J8[. 

In this section first we outline the derivation of this result (for details consult 
refs. (pi ) . Then we generalize it and show that it is correct also for a dense set of maps 
on the boundary of the Mandelbrot set. 

The Julia set J is the closure of all repelling periodic orbits. If c lies outside the 
Mandelbrot set M, the Julia set resembles the Cantor dust. To prove ( |33| ) we have 
to use analyticity properties of the transfer operator and its adjoint. One starts by 
defining the evolution operator £(-i) on the space of functions ip which are locally 
analytic in a small neighbourhood around J . Then one can prove that the spectrum 
of £(—i) is a point spectrum with the only possible condensation at zero. The dual 
or adjoint operator C*,-^ of the transfer operator £(-i), defined on densities 4>* dual 
to the functions ip, then has the same spectrum. The duality is with respect to the 
bilinear form 

W) = / $:<Kx)1>(*) > and 4>{x) = <£*(— ) (34) 

which is natural if one deals with analytic functions. The function <f>(x) is holomorphic 
on the complement of J and vanishes at infinity. The contour 7 encircles the Julia 
set and lies in the common domain of analyticity of <p and ip. The spectral equation 
4>* — zC*,^* — can then be written as 



I ^ {X) 



'7 
for all ip. Therefore 



0(aO - z-±r<Kf(x)) ) = 0. i-r.) 



R(x) := 4>{x) - zj7^j<f>(f( x )) (36) 
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has to be holomorphic around J and can have poles only if ,3 ■, has poles, because 4> 
is analytic outside J. 

Specializing to the map f(x) = x 2 — c, the adjoint equation ( p5| ) leads to 

fl(a;) = 0(a;) - ^(x 2 - c). (37) 

We conclude that R vanishes at infinity because <j> does. Furthermore R can only 
have a pole of order one at zero, the unique critical point x c of /. That means we 
have either R(x) — - (up to a multiplicative constant) or R(x) = 0. R(x) = is 
contradictory because it yields eigenfunctions equal to zero. By iterating ( p7| ) after 
solving it for 4>(x) we get 

i °° n 

<l>(x) = ~ + Y Z —— 1 . (38) 

x r ti2«/™(*)n;Co/^) 

Now note that <j) is only an eigenfunction if it is holomorphic at the poles of R, that 
means at zero in our case (f'(x c ) — 0). Therefore 

= Res x=0 [(f>(x)} 

= 1 + £ 2 n /(0 ) /2( Z 0).../"(0) = F(Z ' ~ X) (39) 



is an equation for the characteristic values of £(-i) meaning that (Bq) is an 
eigenfunction if and only if z is an eigenvalue. For a generalization of this method to 
arbitrary rational maps see S . 

The form ( p3| ) for the spectral determinant for the map f(x) = x 2 — c holds true, 
if c lies outside or inside the Mandelbrot set M ||. However, the method fails if c lies 
on the border of M, and in this case the form of the spectral determinant is unknown. 

As a generalization to || we show here that on a dense subset of the border of M, 
the so called Misiurewicz points, formula ( |33| ) still applies. The Misiurewicz points 
are defined by the requirement that the critical point x c — is preperiodic [O. In 
figure @ we introduce the notation for preperiodic points. If we take the validity of 
( p3| ) for granted, a straightforward calculation yields that for Misiurewicz points the 
spectral determinant is in fact a rational function composed of a transient part and a 
periodic contribution 

* — 2 „ t _i 



F(Z,-1) = 1 + Y * U ., ; + ^ , CT(Z) r ■ (40) 

k ' ' ,ti2«n; = iP(o) 2^n*=lP(o)(i ^ 



A, 



Here 



2 
2xi 4:X\X2 ' ' 2 n p~ 1 xi ■ ■ ■ x r . 



a(z) = l+ + + ... + = , (41) 



A p = 2™ p xi • • • x„ is the stability of the periodic orbit, n p is the period and t the 
transient time of the Misiurewicz point. The n p + t — 2 zeroes of F(z,—1) are the 
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f\0) = xj 




X 3 



Figure 2. The critical point x c = is a preperiodic point with transient time t = 3 
and period n p = 5. In this case the parameter c of the quadratic map is called a 
Misiurewicz point. Note that the periodic orbit is a repeller. 



characteristic values of £(_i)- F has n p poles which lie on a cycle around the origin 
in the complex z-plane. As the period n p of the orbit increases to infinity, these poles 
built up a wall which limits the region of analyticity. 

Presumably, formula ( [f0|) can be proved by a generalization of the arguments in 
p. Here we check it using cycle expansion, i.e. the power series representation (ft) 
of the Fredholm determinant |IJ, [13| . As the first example we compute the spectral 
determinant for the Ulam map f(x) = x 2 — 2. For this map the critical point x c = 
iterates via /(0) = —2 to the unstable fixed point x* = 2. That means we have t = 2 
and n p = 1, and formula (J4CJ) yields 

z 1 1-2/2 
F(z,-l) = l - = V- (42) 

in agreement with the cycle expansion result of ref . E3| , where this formula is derived 
by observing that for all cycles the stability is |A p | = 2™ p , with exception of the fixed 
point xq for which Ao = 4. However, the result for a general Misiurewicz polynomial 
is not as trivial. 

Next we check (J40|) by a numerical computation of the spectral determinant via 
complex cycle expansion. The method is described in detail in |L8| . Here we will only 
state the results. As a first example we study the case c = 1.5436890126920763616 . . ., 
for which the critical point x c = iterates after three steps into a fixed point 
x* = c(l — c), so t = 3 and n p = 1. Therefore we expect a spectral determinant 
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n 


a n 







1.00000000000000000 


3.087378025384153 


1 


-0.3238994356305212 


1.678573510428349 


2 


-0.1929611265865065 


-1.678573510428322 


3 


0.1149554222008833 


-1.678573510428341 


4 


-0.0684839963735331 


-1.678573510428292 


5 


0.0407989259618778 


-1.678573510428285 


6 


-0.0243057129809394 


-1.678573510428318 


7 


0.0144799812638157 


-1.678573510428306 


8 


-0.0086263611178524 





Table 2. Spectral determinant for a chaotic quadratic map up to order 8 in complex 
cycle expansion. For n > 1, the ratio a n /a n +i of the cycle expansion coefficients 
(|) computed from 2™ periodic points equals A = 2c(l - c) = -1.6785735104283 to 
machine precision, in agreement with J43J) 



of the form 



F(z,-l) = l 



2/(0) 4/(0)/2(0) (1 - ^) 



(43) 



c(c — 1). In table g the coefficients a n of the 



n=0 



are depicted. 



with /(0) = -c and / 2 (0) = -x 

cycle expansion for the spectral determinant F(z,— 1) = J^ 

The second column of the table demonstrates that they indeed decrease geometrically 

with the correct ratio. The remaining coefficients can be checked numerically or 

computed analytically from the traces of £(— i) and £>i_-f\- In cither case we arrive at 

(fig). Because the numerator of the spectral determinant is quadratic there are two 

characteristic values, namely z — c 2 * (1 — c) = —1.295597742522084 and z = 2. 

That means the spectrum is real and there is a real pole at z = 2c(l — c) = 

— 1.678573510428321. The escape rate ( pl[ ) is given by the trivial largest eigenvalue 

7 = log(2). 



As our last example we take the map f(x) 
trajectory 



for which x r — has the 



h^ i h^ x\ = —1 + i h^ a; 2 = — i i— > —1 + i i— > . . . , 
so t — 2 and n p = 2. The numerical cycle expansion |L8| yields 



F(z,-1) 



1i 4+4i 



1 



4+4i 



where by (40) we had expected 



(44) 



F(z,-l) = l 



2/(0) 1 



(45) 



4£Ci:E2 
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After short inspection the two expressions turn out to be the same. 

Any number of such examples can be worked out in the same way. In every case 
a scaling of the numerical coefficients confirms the rational form of the Fredholm 
determinant for the Misiurewicz polynomials. 

6. Multi-dimensional polynomial mappings 

In the multi-dimensional case a typical contour integral representation of a periodic 
orbit sum is of form 



/ dx\ / dxi / dx e i -i-i- 1 

J 2m J 2m J 2m ^ f a (x) - x a 



(46) 



For example, in 2 dimensions the contour integral picks up a contribution from each 
fixed point i 



fi(xi, Vi) -Xi=0, h{xi,Vi) - Vi 



0. 



The integral can be converted into a sum of local contour integrals around linearized 
neighbourhoods of the fixed points (xj, yi) 



dzn f dz. 



r(o) = y i 2* a ^ 

^ J 2m J 2m [zn - Jiizn 



1 



Jl2Zi2)(Zi2 — J21ZH — J22Zi2) 

where x = x% + Zn, y = y% + Zi2, and Ji is the Jacobian matrix 



Jn J12 
J 2\ J22 



9xfi(x,y) d y fi{x,y) 
. d x f 2 (x,y) d y f 2 (x,y) _ 



(x,y) = (xi,yi) 



Completing the integrals we get for 2 dimensions 



»eFix(/) 



and in d dimensions 



T(0) 



E 



»ePix(/ 



det(l - Ji) 



)(1 - J22) - J12J21 



0. 



(47) 



where J^ is the Jacobian (monodromy matrix) of the fth fixed point. Pushing all 
contours in ( |46| ) off to infinity yields the sum rule T(0) = 0. Just as in the 1- 
dimcnsional case (2_7|), T(0) vanishes when all fixed points, real and complex, are 
included in the sum. 



6.1. A sum rule related to the fast kinematic dynamo 

Our second multi-dimensional sum rule arises in the study of a 2-d Poincare map of a 
model of the fast kinematic dynamo . The fast kinematic dynamo is a problem of 
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passive vector field advection: the dynamo rate is determined by the overall growth 
of a small vector field embedded into the flow. In this case the transfer operator (g) is 
weighted by the Jacobian matrix J (x, y), a multiplicative function evaluated along the 
trajectory with initial point (x, y). For 2-d maps this transfer operator leads to traces 
of form tr£ = ^ tr J/| det(l — J)|. We can derive a sum rule for the corresponding 
signed measure trace sum: 

T(l)= l dx I dy trJ(ai,y) 



2iri J 2iri (x - f(x,y))(y - g(x,y)) 

tr J, 



E 

iGFix(/ 



det(l 



The 2-d matrix identity 

tr J tr J 



det(l-J) 1-trJ + detJ 

1 1 



det(l-J) det(l-J- 1 ) 



-1 



yields our sum rule 



*W- £ 35^33 -"• (48) 

ieFix(/) v l > 

where we have used the fact that T(0) vanishes for both the forward and the time 
reversed flow. As in the l-d case (Uq), T(l) counts the number of the fixed points of 
the map; in polynomial maps this is given by the order of the polynomial. 



6.2. An application: numerical checks of cycle sums 

As an example for the utility of sum rules for higher dimensional polynomial mappings, 
consider the 2-dimensional Henon map 

x k+1 = 1 - ax\ + y k , y k+1 = bx k ■ (49) 

For the complete repeller case (all binary sequences are realized), the Henon map is 
a realization of the complete Smale horseshoe. Cycle stabilities are easily computed 
numerically [Q . We have verified numerically the sum rules ( f47| ) and ( p8[ ) for a volume 
preserving Henon map repeller (b = — 1) by substituting the cycle stabilities of cycles 
up to length 12 into 

£ X7(i-W = (50) 

^ A^l-l/A,) 2 " ( j 
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(here Aj is the expanding eigenvalue of Jj). The first sum can be used to check the 
accuracy of the periodic orbit data used to compute the escape rate 7 of the Henon 
repeller 

y^ I . ,, -. . . 2 = e~" 7 + (non-leading eigenvalues) . 

ieFix(/™) 

If we add ( |50|) to the above sum, for complete horseshoe all inverse hyperbolic cycle 
contributions with Aj < cancel, and hence only a half of the cycles suffices to to 
compute the escape rate. The second sum rule can be used to check the accuracy of 
the data used to estimate the kinematic dynamo rate r\ 



v A, + 1/A< 

^ |A 4 |(1-1/A0 2 



(non-leading eigenvalues) . 

\ ; I [ I — I / / \ : I ~ 

»eFix(/») 



It was a numerical discovery of these sum rules computed in connection with the 
kinematic dynamo model of ref. |6l that lead us to the contour integral formulation 
and derivation of all other sum rules given in this paper. 

7. Farey map (or spin-chain) thermodynamics 

So far, our examples were based on the contour integral technique. In this section 
we shall use a very different, number theoretic approach to compute the spectrum of 
transfer operators. We shall consider the Farey map, figure H, a combination of two 
Moebius transformations 

Ux)-[ X/{1 ~ X) X G [M/2) (521 

m ~\ (l-x)/x x e [1/2,1] • (52) 

The Farey map arises in a variety of contexts: dynamics of circle maps, dynamical 
renormalization theory, statistical mechanics, number theory. For dynamical systems 
motivation we refer the reader to the review |19| . In particular, the Farey 
map plays the same role vis-a-vis the shift circle map that the highly nontrivial 



circle map presentation function |20|, 21 plays in the golden-mean renormalization 
theory, and provides a controllable setting to test some of the ideas that arise in 
that context. Another motivation comes from statistical mechanics, where Farey 
model thermodynamics corresponds to a number-theoretical model of an infinite 
ferromagnetic spin chain with effective r~ a interaction [ pT[ , with phase transition 
of Thouless type at a — 2. For the purpose at hand the Farey map is of interest to us 
because 

(i) again an infinity of sum rules can be obtained, but this time without the contour 

integration approach employed above 
(ii) the sum rules apply to the positive measure |A| r , in contradistinction to the signed 

measures of the previous examples 
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Figure 3. The two branches of the inverse to the Farey map (p2|). 

(iii) the theory can be extended to no-integer exponents r, with infinite range recursion 
relations. 



Apply the "natural measure" modification of 
£ (T) (y,x) = \f'(x)\ T 5(x-f- 1 (y)) 



(53) 



to the Farey map (we have redefined the exponent (fc + 1) — > r to conform with the 
notation conventions of refs. E3, P3j), and consider the sum 

Z n (r) := f dx \f n '(x)\ T S(x - /~ ( " +1) (l)) . 



(54) 

The leading C eigenvalue 2 9<T ) (in the notation of ref. p^| ) dominates this sum in the 
n —>■ (X) limit and defines a "thermodynamic" function q = q(r). For the Farey map 
(see figure ||) the pre-images of 1 are 

/(1/2) = 1, / 2 (l/3) = / 2 (2/3) = 1 

/ 3 (l/4) = / 3 (2/5) = ■•• = ! 
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f n (Pi/Qi 



i 



These fractions form the levels of the Farey tree |24|, |2^, |2(], |27], |2^, [2£|, a number 
theoretical construction based on the observation that somewhere midway between 
two small denominator fractions (such as 1/2 and 1/3) there is the next smallest 
denominator fraction (such as 2/5), given by the "Farey mediant" (P + P')/(Q + Q') of 
the parent mode-lockings P/Q and P'/Q'. The Farey tree is obtained by starting with 
the ends of the unit interval written as 0/1 and 1/1, and then recursively bisecting 
intervals by means of Farey mediants. This generates Farey level sets T n with 2 n 
mode- locking widths on each level, figure 0. In the present context, the nth level of 
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Figure 4. The levels T n of the Farey tree are generated by the Farey mediant 
addition rule P/Q, P'/Q' -» (P + P')/(Q + Q')- 



the Farey tree is the set of the 2 n distinct backward iterates f~ n (xi) = 1 of the Farey 
map. Furthermore, noting that f' (x) — 1/(1 — a;) 2 , f[(x) = —1/x 2 , it is easily checked 
that if Pi/Qi € T n , then 



\f n '(P i /Qi)\=Q 2 i 



(55) 



For example 



/ 2 '(l/3) = /'(l/3)/'(l/2) = (1 - l/3)- 2 (l - 1/2)- 2 = 3 2 



Hence the sum (B4) is the sum over denominators of the Farey rationals of the (n-l)th 
Farey level: 



Z n {r) = ]T QT- 



(56) 



ier„ 



where Qi is the denominator of the «th Farey rational Pi/Qi- For example 
Z 2 (l/2) = 4 + 5 + 5 + 4. 
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As we shall now show, Z n (r) satisfies an exact sum rule for every non- negative 
integer It. First one observes that Z n (0) = 2". It is also easy to check that J24[ 
Z„(l/2) = YsiQi = 2 - 3". More surprisingly, Z„(3/2) - £ 2 Q 3 = 54 • 7"" 1 . Such 
sum rules are consequence of the fact that the denominators on a given level are Farey 
sums of the denominators on preceding levels. In order to exploit this, the following 
labelling of the Farey denominators introduced by Knauf p9| is convenient: 

The Farey denominators Q a , a <E {s\S2 ■ ■ ■ s n }, Si € {0, 1}, up to the nth level 
of the Farey tree of figure [| can be labelled in lexical binary order as illustrated in 
figure 0. Farey denominators on the (n + l)th Farey tree level are given by p9[ 



n = 




Figure 5. Lexically ordered binary labelling of all Farey denominators up to and 
including the nth Farey tree level, following Knauf |2S|]. The binary label of the 
denominator Q CT is read off starting with the root of the tree, with Sj = 1 for a full 
line, and Sj = for a dotted line. For example Qooil = 7, Qiooo = 2, etc.. The 2 n 
denominators Q.i preceded by a full line constitute the nth Farey level T n - 



QoQ = Qa (carry over all previous Farey denominators) 

Qui = Qa + Qw (compute the (n + l)th Farey tree level) , (57) 

where a is obtained by exchanging l's and 0's in the string a, Si — ► 1 — Sj, and 
the recursion is initiated with the empty string values Q, = Q_ = 1 . The second 
rule follows from labelling of nearest neighbours on a given level, see figure |5j. By 
construction, Q a \ is invariant under the s, — > 1 — s, interchange 



Qcrl — Q~5\ 



(58) 



Define a [2r + l]-dimensional vector <f> = [<j>o, fa, ■ ■ ■ , far], It a positive integer, by 
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The zeroth component of <p^ is the /cth level Farey sum (|5 
Z k (r) = ^= E <&■ 

o'G{siS2---Sfc} 



(59) 



Motivation for constructing the vector (j^ k > is the observation that its remaining entries 
exhaust all combinations that arise in exponentiating Q a \ = Q a + Q-& to the integer 
powers s < 2r 

<r£{siS2---s k+1 } 
r?£{sis 2 ---s k } 

= E w 

o-e{sis 2 ---s fe } 



/r. ] 1- Qcro) r m (Q™1 + Q™l) 



(we have used ^ CT = J^— and the Qo-i = Qwi symmetry), and that they form a linear 
basis for d>^ k+1 ^>: 



2r—m 

€ +1) = E 

r=0 

It— m 



E 

r=0 



2r — m 
r 

2r — rn 
r 



<TG{S1S2---S fe } 



2t — TO 



4 fe) . 



Hence the family of vectors <p l \ fi^K ' ' '> is generated by multiplication <j)( k+1 
L( r )0( fc ' by the [(2t + 1) x (2r + 1)] dimensional transfer matrix 

/2t — to\ /2t — to 

L(2r)mr = + 

\ r y \ r — to 
consisting of two Pascal triangles. For example 
2 4 2 



(60) 



L (2) - 



1 2 1 
1 1 



The growth of (f)^ k > in the k — > oo limit (and in particular, its zeroth component (|59|), 
the thermodynamic sum Zk(r)) is given by 1/z = 2 <? ( T - ) , the leading zero of the L(2r) 
characteristic polynomial 

F{z, 2r) := dct (l - zL (2t) ) = 1 - (2r + 1 + F 2r+1 )z - . . . , (61) 

see table g. Here F n = F„_i + F n -%, Fq = 0, F\ = 1, are the Fibonacci numbers. The 
largest Qi in the nth level sum (|56|) is -Fn+2, so for the large r the leading eigenvalue 
tends to p 2r . The polynomials of table [| were first computed by Cvitanovic and 
Kennedy, Appendix B of ref. 1 23 ; here we have followed the more elegant derivation 
due to Contucci and Knauf flll|. 
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Table 3. Spectral determinants for the Farey model transfer operator 

It = 0, 1, 2, ... , 11, together with the leading eigenvalue 2 q ( T >. 



2r 



2«( r ) 



F(z,r) 



1- 2z 
l-3z 

l-5z + 2z 2 

l-7z 

(1 + z) (l- llz + 2z 2 ) 

1- 14z-47z 2 

1 - 20 z - 161 z 2 - 40 z 3 + 4 z 4 

1- 29z-485z 2 -327z 3 

(1 + z) (l - 44 z - 1313 z 2 - 88 z 3 + 4 z 4 ) 

1 - 65 z - 3653 z 2 - 3843 z 3 

(1 + z) (1-z) (l - 100 z - 9601 z 2 - 200 z 3 + 4 z 4 ) 

1 - 156 z - 24882 z 2 + 83828 z 3 + 107529 z 4 

(1 + z) (l - 247 z + 63659 z 2 + 797003 z 3 - 127318 z 4 - 988 z 5 + 8 z 6 ) 
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(1.62413- 


.)22 



in the — > oo limit, 



(l + V5)/2 =1.61803. 



The Farey model demonstrates that not only are the methods of this paper 
competitive with the periodic orbits methods, but they are sometimes superior to 
them. For the Farey model the periodic orbit expansions are studied in detail in 
refs. [ p0[ pl|, and depend on quadratic irrationals rather than the Farey denominators. 
While such periodic orbit expansions are analytically intractable, this variant of the 
finite level sums yields q(r) exactly for all 2r a non-negative integer. Not only that, 
but Contucci and Knauf JTTJ have in this case been also able to analytically continue 
the L( T ) matrices to arbitrary real positive r. 

8. Some speculations 

So far we have derived a series of exact sum rules for very specific simple dynamical 
systems. In this section we speculate about possible, but at present still largely 
unexplored applications of the above methods in more general settings. 



8.1. Normal form approximations to flows 

The Poincare maps of flows of physical interest are in general not polynomials, but 
smooth analytic functions. However, in practical applications (such as the long term 
integrations in celestial mechanics) it is often advantageous to replace a differential 
equation by a normal form that approximates the return map for a Poincare section 
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of the flow, and replace the numerically demanding integration of the flow by a 
map iteration. Such applications require a test of the accuracy of the normal form 
approximation; the above sum rules offer an estimate of the quality of approximation, 
which, as it depends only on cycle eigenvalues, is coordinatization independent. The 
idea is simple: were the polynomial representation of the flow exact, its cycles would 
satisfy the above sum rules exactly. By comparing the sums of cycles of the smooth 
flow to what would be expected were the flow generated by a polynomial we learn how 
good a polynomial approximation of given order would be, prior to any actual fitting 
of the flow by an approximate mapping. 

As an example, consider a 3-dimensional flow x = F(x), x = (xi,X2,x 3 ). We 
assume that the flow of interest is recurrent, and that given a convenient Poincare 
section coordinatized by coordinate pair [x,y), the flow can be described by a 2- 
dimensional Poincare map 

y' = g(x,y) 

together with the "ceiling" function T(x, y) which gives the time of flight to the next 
section for a trajectory starting at (x,y). In general / and g can be complicated 
functions, but the essential properties of a continuous flow can be modelled by the 
Henon map (|49|), which we take as a local normal form (up to quadratic terms). In 
the Hamiltonian case, the normal form is 

x k+1 = 1 - ax\ - x fe _i . (63) 

If this were an exact representation of the flow, it would satisfy exactly the sum 
rules ( pOJ ) and (pl|), hence the amount by which they are violated is an indication of 
how well the flow is approximated by a quadratic normal form. As we are summing 
signed rather than absolute values, this is an underestimate of the actual error, but 
possibly of interest even so. While numerical tests of the rule on a 3-disk billiard cycles 
data of ref. JL2) do suggest that for sufficiently large disk-disk separation Jp3j ) is not a 
bad approximation, we have not yet checked whether the best value of the parameters 
in the approximation can be determined by minimizing the deviation from the exact 
sum rules. 

8.2. Spectrum of the period- doubling operator 

The generalized transfer operators (||) have an important application |IJ, |3^, [33| Q 
in the period-doubling renormalization theory: if / is given by the period-doubling 
presentation function | p8[ 

f (x) = ag(x), g (a" 1 ) < x < I 
fl(x) — ax, aT x < x < a~ 2 , 

where 

g(x) =agog(x/a) . 
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is the universal period-doubling renormalization fixed point function, the leading 
eigenvalue of C(2)(y,x) = S(y — f^ 1 (x))f'(y) is the Feigenbaum S, and the spectrum 
of £(2) is given by the trace formula pi p3| 

tr£ ( 2 )= E a^T' (64) 

i6Fix(/) l 

or the associated Fredholm determinant. Here the summation is restricted to the 
real fixed points in the unit interval; application of our sum rule (|19| ) would require 
inclusion of the complex periodic points as well. The crudest approximation to g{x) 
is a quadratic polynomial. This already yields a reasonable estimate without any 
actual polynomial fitting; according to (191), for a quadratic polynomial the leading 
eigenvalue (pit) is 4, while the period-doubling leading eigenvalue is the Feigenbaum 
8 = 4.66 .... If (|19| ) is really applicable in this context, the sum rule would require 
inclusion of contributions of complex periodic points of the presentation function - 
how well that converges has not yet been checked. 

9. A critical summary 

The conventional periodic orbit theory suffers from one serious limitation; the number 
of spectral eigenvalues computable is in practice limited by the exponential growth 
in the number of periodic orbits required. This arises because the periodic orbit 
theory in its trace formula formulation is essentially a tessellation of the phase space 
by linearized neighbourhoods of periodic points which ignores the fact that these 
neighbourhoods are interrelated by the analyticity of the flow. This suggests that 
the local periodic orbit information should be supplemented by global analyticity 
constraints. Previous to the work presented here, analyticity has been exploited to 
improve the convergence of cycle expansions by incorporating the shadowing of long 
cycles by combinations of shorter ones J35| [l4| , and to prove that the cycle expansions 
of Fredholm determinants can converge faster than exponentially with the topological 
cycle truncation length |3q, p7|. While sometimes very useful, these results apply 
only to nice hyperbolic flows with finite symbolic dynamics. Most realistic physical 
flows arc not uniformly hyperbolic, and do not have finite Markov partitions. For 
this reason general analyticity constraints applicable to any flow would be of great 
interest. Here another observation is very suggestive; analytic continuation of pairs of 
orbits through inverse bifurcations J38| improves the semi-classical quantization in a 
regime where the dynamics is non-hyperbolic, and utility of periodic orbits a priori in 
doubt. Similarly, leading diffraction effects of wave mechanics can be accounted for by 
inclusion of complex "creeping" orbits |39| ] into periodic orbits sums. These examples 
suggest that for generic flows perhaps both the non-hyperbolicity and lack of finite 
Markov partitions may be tamed by taking into account the totality of periodic orbits, 
real and complex. 
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With such considerations in mind we have here compared two strategies to compute 
the spectrum of a given transfer operator: 

(i) the conventional periodic orbit theory approach: tr C = sum over cycles 
(ii) the new recursive approach: relate Z n = J dx £ n (0, x) = sum over pre- images to 
previous levels sums. 

We have established that sometimes the iterative nature of the dynamics enables us 
to obtain formulas for transfer operator eigenvalues without computing any periodic 
orbits. 

The sum rules so obtained come as something of a pleasant surprise, especially 
in the multi-dimensional cases, as cycle expansions typically require evaluation of 
exponentially many cycles up to a given cycle length, and it is not at all obvious 
that the cycle weights should be related in any simple way. Prior to these results, 
other than the purely topological Lefschetz fixed point formulas there existed only 
one exact general result for periodic orbit sums, the flow density conservation J4fJ Q, 
and several sum-rules specific to billiards. These are explored in a forthcoming paper 
by S.F.Nielsen et al. M). From the vantage point of this work, the flow density 
conservation rule is of limited interest, as it has no dynamical content. Even though 
our sum rules require determining not only real cycles, but also the complex ones 
(ie. also those forbidden by the symbolic dynamics), they might be useful in practice. 
There are several immediate uses for such rules; one is to check for the correctness 
of the numerically calculated cycle weights up to a given cycle length, and another, 



speculated on in section 8.1, is to estimate the quality of polynomial approximations 
to smooth flows. 

Applied to higher dimensional dynamical systems the new sum rules might shed 
light on the role of complex periodic orbits in phenomena such as intermittency 
and pruning, reduce the number of periodic orbits required for accurate eigenspectra 
computations, provide useful cross checks of the accuracy of trace formulas and cycle 
expansions, might lead to replacement of periodic orbit sums by contour integrals, and 
provide error bounds on approximations of smooth flows by normal form mappings. 
Still, the explicit sum rules for simple models presented above are only a modest step 
in direction of moving "beyond periodic orbit theory" : 

(a) The results presented depend on either the finiteness of polynomials (for the 
contour integral method) or properties of integers (for the Farey map case). One 
needs to develop infinite dimensional versions of the above finite rank operators 
L (ref. Q is a step in this direction). 

(b) One needs to reinstate positive measures. While in quantum mechanical 
applications the signed measures might be just the right Maslov phases, in the 
classical context the above sum rules are essentially index theorems and play a 
role analogous to what Lefschetz zeta function plays vis-a-vis the Artin-Mazur 
counting zeta function. (Appendix of ref. |9| is a step in this direction) . 
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(c) The transfer operator weights studied here were restricted to powers of the cycle 
stability; one needs to deal with more general classes of weighted operators. 
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Appendix A. Examples of sum rules 

To give the reader a feeling for what the above sum rules mean in practice, we shall 
work out explicitly the few simplest examples. The quadratic polynomial fixed points 
and their stabilities are given by 



x = (l + Vl + 4c)/2, Xx = (1 - VI + 4c)/2 , 

A = 2x , A x = 2xi , (Al) 

and the 2-cycle stability by 
A i = 4(1 - c) . 

Appendix A.l. Trace formula for Ci \ 

The sum rule dlj) states that for the quadratic polynomial 



11 11 

= - - + - -, = 



Ao-1 Ai-1' Ag-1 Af-1 Aoi-1' 

and so on. This sum rule is not new; it goes back to Fatou G3 and Julia [E5[ and can be 
found, for example, in sect. 3.5 of ref. |15| , proven by manipulating roots of polynomials 
rather than by the contour integration method. The Fredholm determinant (|2(]) is in 
this case trivial, F(z, 0) = 1. 
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Appendix A. 2. Trace formula for C(y\ 

The sum rule (Uq) for the quadratic polynomial fixed points and 2-cycles can be 
checked by hand: 

o_ A , A l „ _ A , A ? ,. 2 _ A 01 



Ao-1 ' Ai-1' Ag-1 A?-l Aoi-1 

For quadratic polynomials TV = 2™, so both the Fredholm determinant (gQ) and the 
dynamical zeta function (pG) are given by F(z, 1) = 1 — 2z, 1/C(i) = F(z, l)/F(z, 0) = 
1 — 2z . For the case of all periodic points real, the traces T n (l) count the numbers of 
periodic points of period n, and 1/C(i) is the Artin-Mazur topological zeta function Eg . 

Appendix A. 3. Trace formula for C(2) 

The sum rule (|19|) for the quadratic polynomial is given by 

A 2 

^ A-l 
teFix(/») l 

the Fredholm determinant by F(z, 2)(z) = 1 — 42, and the dynamical zeta function by 

F(z,2) l-4z 



VC(2) 



F(z,l) l-2z 



Furthermore, as for every k > the term rflQ ) contributes a factor (1 — 2 fe z) to F(z, fe), 
1/C(fc) always contains a factor (1 — 2 fe z)/(l — 2 k ~ 1 z), in agreement with ref. [jHj. 

Appendix A. 4- Trace formula for £( 3 ) 

According to (|1J), we only need to to evaluate C(3) for the quadratic polynomial 

^,„s f dx 8x 3 „ f dx ,„ c , 

(7 3 = ^^-^ = 8 0— ar(l + -j + ...)= 8c. 

J 27TI £ — C J Z7TI X z 

Together with (pin this yields the T„(3) sum rule for the quadratic polynomial 

E a^~i = 8 " + (8c) "- (A2) 

ieFix(/») ' 
For example, it is easily verified that 

Ag A? 



1 + c = 



A| A^j A| 

a§-i + a 2 -i + 'a 01 - 



Ao - 1 Ai - 1 

64 < 1 + c2 ^aF L T + a& + 2 

The Fredholm determinant is given by 

F(z,3)(z) = (l-8z)(l-8cz) , 
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in agreement with the k = 3 case of our general formula (|27J). The dynamical zeta 
function is 

1; , F(z,3) (l-8z)(l-8cs) 

1/C ( 3 > = fM = l^iz ' (A3) 

Thus the first "nontrivial" zeta function is indeed of the form conjectured in ref. p0[ . 

Appendix A. 5. Trace formula for C(j\ 

We finish this appendix by writing out explicitely a typical nontrivial example: 
k = 7 = 2x3 + 1. According to (||) 

■c 3 3c 



L(7) 



c 4 6c 2 1 
c 5 10 c 3 5 c. 



The corresponding Fredholm determinant is given in table El This, as well as any 
other F(z, k) that we have computed agrees with the form conjectured in ref. ]10[ . 
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